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ABSTRACT 

Using a set of compilations of measurements for extragalactic radio sources we con- 
struct all-sky maps of the Faraday Rotation produced by the Galactic magnetic field. 
In order to generate the maps we treat the radio source positions as a kind of " mask" 
and construct combinations of spherical harmonic modes that are orthogonal on the 
masked sky. As long as relatively small multipoles are used the resulting maps are quite 
stable to changes in selection criteria for the sources, and show clearly the structure of 
the local Galactic magnetic field. We also suggest the use of these maps as templates 
for CMB foreground analysis, illustrating the idea with a cross-correlation analysis 
between the Wilkinson Microwave Anisotropy Probe (WMAP) data and our maps. 
We find a significant cross-correlation, indicating the presence of significant residual 
contamination. 
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1 INTRODUCTION 

The origin of large-scale magnetic fields observed on galac- 
tic and cluster scales is unknown. The magnetic fields, with 
observed strengths of ~ 10 _6 G, could be the consequence 
of an amplification of a tiny seed (5s 10 _20 G) by a dynamo 
mechanism. Alternatively, the compression of a primordial 
seed (~ 10 _9 G) by protogalactic collapse could lead to the 
fields we see today. Both scenarios require an initial pri- 
mordial field. Furthermore, the two mechanisms need to ex- 
plain the high redshift magnetic fields observed in galaxies 
(Kronberg, Perry & Zukowski 1992) and damped Lyman-a 
clouds (Wolfe, Lanzetta & Oren 1992). Magnetic fields play 
a crucial role in star formation as well as possibly playing 
an active role in galaxy formation as a whole (Wasserman 
1978; Widrow 2002). Thus, one of the most significant tasks 
in cosmology is unravelling the mystery behind magnetic 
fields; from the primordial field to our own Galactic field. 

The cosmic microwave background (CMB) provides us 
with the most distant and extensive probe of the early uni- 
verse. A primordial magnetic field will leave an imprint in 
this radiation. Various methods have been developed that 
seek these signatures. Barrow, Ferreira and Silk (1997) use 
the anisotropic expansion caused by the presence of a ho- 
mogeneous primordial field to place limits on its size from 
large- angle CMB measurements. Others have sought correla- 
tion between different scales in the temperature anisotropies 
(Chen et al. 2004; Naselsky et al. 2004) or computed the 
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effects the field has on the polarisation-temperature cross- 
correlation (Scannapieco & Ferreira 1997; Lewis 2004). The 
existence of a magnetic field at last-scattering also leads to a 
possible measurable Faraday rotation of the polarised CMB 
light (Kosowsky & Loeb 1996). 

At the other end of the scale, investigation of our own 
Galaxy's magnetic field has led to the development of a num- 
ber of different techniques; Han (2004) provides a concise 
review of the subject. Techniques include observing Zeeman 
splitting, polarised starlight, synchrotron radiation, Faraday 
rotated light and polarised dust emission. However, even 
with all these methods, there are still outstanding problems. 
This has led to a lack of consensus on key issues: from the 
number of spiral arms; how the arms are connected; to the 
direction the field takes along the arms (Vallee 1997; Han 
2004). 

To fully understand magnetic fields, we need a coherent 
picture throughout different epochs. Theories and models 
can then be tested against this observational picture. A ro- 
tation measure (RM) map of the full sky has the potential 
to fulfil such a goal. RM values probe the integral of the 
magnetic field from the radiation source to the observer. 
Obviously, the information encoded in such a map depends 
on the location of the radiation that is rotated. Faraday ro- 
tated polarised CMB radiation will offer both a picture of 
the primordial field (at the surface of last scattering) and 
that of our Galaxy. The recent detection of CMB polarisa- 
tion by DASI (Kovac et al. 2002) and confirmation of this 
via the Wilkinson Microwave Anisotropy Probe (WMAP; 
Kogut et al. 2003), have opened up a new avenue in CMB 
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research. Future results from WMAP and Planck satellites 
will offer polarised data covering the full-sky over a range 
of frequencies (seven in the case of Planck). By forming a 
RM map with the data and looking at differing scales, we 
should be able to untangle the information in the data; on 
the largest scales the local magnetic field can be studied, 
whereas the primordial field can be studied on the smaller 
scales. 

On the other hand, one of the main tools for probing 
local magnetic fields (such as our Galaxy's) involves utilising 
RM values from extragalactic sources. Catalogues contain- 
ing RM values of extragalactic sources have been used to 
map the Galactic field (eg. Frick et al. 2001). Combining this 
data with rotation measures from pulsars that are located 
within our Galaxy, it is conceivable that a 3-dimensional 
image of the Galactic magnetic field can be built. However, 
current RM catalogues are both sparsely populated and un- 
evenly sampled. Thus, astute methods are required to pro- 
duce a RM map with the data at hand. 

So, what do we intend to do? We attempt to map the 
RM values as a function of angular position giving 1Z(Q,), 
where we use 1Z to denote the Faraday rotation measure and 
Q for the angular position. Catalogues containing RM values 
of extragalactic sources are used to construct the function. 
The observed spatial distribution of the RM values can be 
expanded over a set of orthogonal basis functions. For anal- 
ysis of data distributed on the sky, expansion over spherical 
harmonics seems natural 

oc m— + Z 

TZ{Q) = J2Yl a '^Yi, m (n), (1) 

Z— 1 m— — I 

where the ai <rn are the spherical harmonic coefficients and 
the Yj jm are the spherical harmonics. The properties of the 
spherical harmonics are well understood and the calcula- 
tion of the ai }m will allow us to utilise routines within the 
HEALPixt package (Gorski, Hivon & Wandelt 1998) for vi- 
sualisation purposes and further analysis. However, a non- 
uniform distribution of data points compromises spherical 
harmonic analysis due to the loss of orthogonality (Gorski 
1994). It is more fruitful to analyse a system using orthogo- 
nal functions: the statistical properties of the coefficients are 
simplified. If nonorthogonal functions are used, the proper- 
ties of the system and the basis are confused. Therefore, 
we would like to construct an orthonormal basis with func- 
tions closely related to those of the spherical harmonics. The 
spherical harmonic coefficients can then be obtained from 
the resultant coefficients of the orthonormal basis. 

Spherical harmonic analysis of extragalactic sources has 
been previously performed by Seymour (1966,1984). How- 
ever, the analysis was carried out using a different form of 
orthogonalisation and only on a set of 65 sources. 

The RM map resulting from our method will be a use- 
ful tool for probing Galactic magnetic structure. The map 
will also be a valuable point of reference when investigating 
mechanisms that involve the Galactic magnetic field. For ex- 
ample, CMB foregrounds (synchrotron, dust, free-free emis- 
sion) are correlated with rotation measures (Dineen & Coles 
2004). There is evidence of another foreground component 
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(de Oliveira-Costa et al. 2004), labelled foreground X, that 
is spatially correlated with 100/im dust emission. Spinning 
dust grains (Draine & Lazarian 1998) are the most popu- 
lar candidates for causing this anomalous emission. A RM 
map can provide insight into the role of these grains as they 
may align with the local magnetic field. Foreground removal 
will be particularly challenging in CMB polarisation stud- 
ies as foregrounds are more dominant than in the temper- 
ature anisotropies. Also, single frequency polarisation mea- 
surements will not be able to remove the effects of Faraday 
rotation through the Galactic magnetic field. Thus, the ex- 
tent to which the results have been effected by the E-mode 
signal rotating into the B-mode signal (and vice versa) is 
unknown. 

The layout of the paper is as follows. In the next sec- 
tion we describe the three rotation measure catalogues used 
in our analysis. In describing the data we clarify the mean- 
ing of extragalactic rotation measures. In Section 3 we il- 
lustrate a method to generate orthonormal basis functions 
for each catalogue. From the coefficients of the new basis, 
the spherical harmonic coefficients are calculated. In Sec- 
tion 4 we present the resulting RM maps and discuss the 
observed features. In Section 5 we give a brief application 
of the maps. Correlations are sought between the RM maps 
and cleaned CMB-only maps. The conclusions are presented 
and discussed in Section 6. 



2 ROTATION MEASURE CATALOGUES 

Faraday rotation measures of extragalactic radio sources are 
direct tracers of the Galactic magnetic field. When plane- 
polarised radiation propagates through a plasma with a com- 
ponent of the magnetic field parallel to the direction of prop- 
agation, the plane of polarisation rotates through an angle 
4> given by 

4> = n\ 2 , (2) 

where the Faraday rotation measure is measured in radm~ 2 
where 

Note that B\\ is the component of the magnetic field along 
the line-of-sight direction. The observed RM of extragalac- 
tic sources is a linear sum of three components: the intrinsic 
RM of the source (often small); the value due to the in- 
tergalactic medium (usually negligible); and the RM from 
the interstellar medium of our Galaxy (Broten et al. 1988). 
The latter component is usually assumed to form the main 
contribution to the integral. If this is true, studies of the 
distribution and strength of RM values can be used to map 
the Galactic magnetic field (Vallee & Kronberg 1975). Even 
if the intrinsic contribution were not small, it could be ig- 
nored if the magnetic fields in different radio sources were 
uncorrelated and therefore simply add noise to any measure 
of the Galactic field (Frick et al. 2001). In a similar vein, 
the distributions of RM values have been used to measure 
local distortions of the magnetic field, such as loops and 
filaments, and attempts have also been made to determine 
the strength of intracluster magnetic fields (Kim, Tribble & 
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Kronberg 1991). In what follows we shall use RM values ob- 
tained from three catalogues in an attempt to map TZ(Q) 
over the whole sky. 

All three catalogues are sparsely populated and have 
non-uniform distributions. It is essential to remove the struc- 
ture due to the spatial distribution of the sources. This 
structure is unique to each catalogue. Therefore, for each 
catalogue a new set of orthonormal functions has to be gen- 
erated. 

The three catalogues used are those of Simard- 
Normandin et al. (1981; hereafter S81), Broten et al. (1988- 
updated in 1991; B88) and Frick et al. (2001; F01). S81 
present an all-sky catalogue of rotation measures for 555 
extragalactic radio sources (ie. galaxies and quasars). B88 
and F01 contain 674 and 800 sources respectively. In F01, 
the two other catalogues are combined with smaller studies 
of specific regions in the sky (see paper for details). They 
also provide slightly reduced versions of the other two cat- 
alogues. Sources with significantly larger RM values than 
those in the other studies are removed, leaving catalogues of 
551 sources for S81 and 663 for B88. In our analysis we will 
use these versions of the two catalogues. 

Finally, we reject sources with TZ > 300radm~ 2 . Such 
large RM values are unlikely to represent real features of the 
Galactic magnetic field: probably they arise from incorrect 
determination of TZ due to the nn ambiguity in polarisation 
angle; magnetic fields within the sources; Equation (2) being 
incorrect; and so on. This final selection criteria reduces the 
catalogues to 540 sources for S81, 644 for B88, and 744 for 
F01. 



3 GENERATING A NEW BASIS 

We can only observe RM values where there happens to be 
a line of sight. This means we see the RM sky through a pe- 
culiar "mask" . We wish to generate a new orthogonal basis 
that takes account of the spatial structure of this mask. In 
particular, we need to find a set of functions that are orthog- 
onal on the incomplete sphere. Ideally, these new functions 
should be related to the spherical harmonics (which are or- 
thogonal functions on a complete sphere). This will enable 
us to determine the spherical harmonic coefficients from the 
new functions and their coefficients. 

Gorski (1994) tackles the problem from the point of 
view of CMB analysis. The determination of the angular 
power spectrum is a crucial element of much work in the 
field. If the temperature anisotropics form a Gaussian ran- 
dom field then they can be completely characterised by the 
angular power spectrum. In order to obtain the angular 
power spectrum, one needs to estimate the spherical har- 
monic coefficients. At low Galactic latitudes (b < 20°) fore- 
ground contamination is severe. Therefore, it is preferable to 
obtain an estimate of the spherical harmonic coefficients out- 
side this region. Gorski (1994) calculates a new set of func- 
tions that are orthogonal to this cut sphere. These functions 
are used to calculate the spherical harmonic coefficients and 
thus estimate the angular power spectrum from the two-year 
COBE-DMR data (Bennett et al. 1994). 

In order to see how the method works it is prudent to 
look at the definition of orthogonal functions. Let us consider 
two complex functions A(x) and B(x). If 



/ 

J a 



A*(x)B(x)dx = 0, 



(4) 



then A(x) and B(x) are orthogonal over on the interval 
{a, b}. If we incorporate these two functions into a vector 
v=[A(x),B(x)], the orthogonality of the functions can be 
expressed through the scalar product 



(v- v T ) {a , 6} 



(5) 



We shall now look at orthogonal functions in the con- 
text of a complete sphere. A function can be de- 
scribed by spherical harmonics Yj jm (0) up to an order 
/max- We can form an (Z max + l) 2 -dimensional vector 

y=[F ,o(fi), Ki,-i(fi), Yi,o(fi), ■ ■ ■ . *i,»«,im«(fi)]- 

The scalar product is then defined as 



(y • y T ) {mi sky> = i- 

The sky can therefore be fully described by 



*+ir 



TZ(Q) = a i Y i(ty =a T y. 



(6) 



(7) 



However, when the sphere is incomplete due to a Galaxy cut 
or a more complex mask being applied, we have 



(y ■ y T ){c UtS ky} = w / 1, 



(8) 



where W is the coupling matrix. A new basis, where the 
equality is true, can be constructed in the following manner. 
The procedure is a type of Gram-Schmidt orthogonalisation. 
W can be Choleski-decomposed into a product of a lower 
triangular matrix L and its transpose 

W = L L T . (9) 

The inverse matrix T — L 1 is then computed. The new set 
of functions on the cut sky is 



1p — (r • y){cut sky}- 

By construction, we have 

(V> • V' T ){cutsky} = r-y-y T -r T 
= r l l t r T 



(10) 



= L 
= I. 



1 • L • L T • (L _1 ) T 



(11) 



This can be a useful cross-check for testing the code. Finally, 
the new basis functions can be used to describe TZ 



TZ(Q) 



J2 



•V>- 



(12) 



At this point, only the angular positions Q n of the 
sources have been required. In order to obtain the coeffi- 
cients a of the new basis, the rotation measure TZ n them- 
selves are required. It is worth putting it into the context 
of the data we have. Let the number of sources in our cat- 
alogue be N and let (! maj + l) 2 =M. It should be evident 
that we have a set of simultaneous equations which have to 
be solved in order to obtain the coefficients of the new basis 
functions 

ft(fil) = Clll>l(Sll)+C2ll>2(Sli) + ...+CMlpM(tol) 

TZ(Q.2) = cn/)i(n 2 ) + 02^2(^2) + • ■ ■ + cmV'm^) 
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TZ(fl N ) = ci^i(fljv) + c 2 ^2(^n) +-.- + CM^M(Hiv).(13) 

So, as long as TV > M, these equations should be solvable. 
Ultimately, we wish to obtain the spherical harmonic coeffi- 
cients ffl;, m . Using Equation (10), we see 

TZ = a T ■ y = c T ■ V = c T ■ T ■ y (14) 

and therefore 

a T = c T ■ T -v a = T T c. (15) 

So, we have obtained the spherical harmonic coefficients. 

There are some practical points that have been glossed 
over in the above description of the method. Firstly, the RM 
values in the catalogues need to be smoothed. Otherwise, as 
the series in Equation (12) is finite, we will be attempting 
to fit large-scale waves to small-scale features. Ideally, the 
smoothing will take place in the new basis, however, this 
is impractical. Therefore, we chose to smooth in harmonic 
space. Around each source, a hoop of 20° is thrown and the 
average RM value is taken of the sources within the hoop. 
This could have been done via a more sophisticated method, 
say a Gaussian-weighted mean of RM values. However, we 
chose to use the simple approach. The size of the hoop was 
chosen to match £ max (~ 16) closely in angular size. Further- 
more, it coincided with the limiting resolution of the wavelet 
method used in Frick et al. (2001) on the same catalogues. 

Secondly, we need to determine to what order we take 
the series up to, i.e. the value of Z max - We do this through 
trial and error. As we will expand upon in the next section, 
RM maps were generated for Z max values of 8, 10, 15, 16, 
17 and 18. The power spectrum for each map was studied. 
At some limiting value of Z max the shape at low I alters 
as features become unstable. Maxima and minima points 
explode since we are trying to fit more and more function 
to the same amount of data. That is to say, N is getting too 
similar to M. 

Finally, the convention for spherical harmonics has to 
be chosen carefully. The new basis functions were calculated 
using the convention of Gorski (1994) where 



O 




Angular Scale I 

Figure 1. Power spectra (measured in rad 2 m~ 4 ). Red i max =15, 
green Z max =16, blue Z max =17, dotted Z max =8 and 10, and dash- 
dot Z max =18. From top to bottom: S81, B88 and F01 catalogues. 



and f(cf>) = tt- 1/2 cos(m0), (2tt)- 1/2 , or tt~ 1/2 sin(|ra|0) for 
m > 0, = 0, or < 0. The spherical harmonic coefficients 
can then be trivially converted into those that adhere to the 
HEALPix definition of the harmonics where now /(</>) = 
(27r)- 1/2 [cos(m0) + isin(m0)] or (-l) m (27r)* 1/2 [cos(m0) + 
i sin(m</))] for m > or < 0. The convention of Gorski (1994) 
was chosen since the information within the coefficients is 
more highly compressed. Whereas in the HEALPix defini- 
tion, the coefficient are complex with a symmetry between 
+m and -m, those following the convention of Gorski (1994) 
are real and contain no such symmetry. The redundant in- 
formation in the HEALPix coefficients leads to confusion 
when solving Equation (13). 



4 ALL-SKY RM MAPS AND THEIR 
INTERPRETATION 

For all three catalogues, sets of spherical harmonic coeffi- 
cients were calculated with Z max being set to 8, 10, 15, 16, 
17 and 18. From these coefficients, RM maps were produced 
using the 'synf ast' routine in the HEALPix package. To see 
whether a RM map was displaying real features or whether 
the series expansion had been extended too far, we looked 
at the angular power spectrum of each map. The angular 
power spectrum is the harmonic space equivalent of the au- 
tocovariance function in real space. It is defined as 

c < = 2TTl£l a -l 2 - (17) 

m 

The spectra of the RM maps are shown in Figure 1. For all 
three catalogues, it is clear that the shape of the spectra is 
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Figure 2. S81 catalogue with 540 sources. Top: Z m ax =10. Bot- 
tom: i max =16. All maps arc shown in Galactic coordinates 
with the Galactic centre in the middle and longitude increasing 
from right to left. The temperature-colour scales arc measured in 
rad m~ 2 . 



consistent up to Z max = 16. Extending the series expansion 
to higher values of / max leads to fluctuations in the power 
on the largest-scales (low I). The method finds it harder 
to reconcile the data with the increasing number of basis 
functions. The result is that maxima and minima explode 
as too much freedom is given. This is clearly visible in the 
maps for Z m ax = 17 and 18 (not displayed). Although looking 
at the bottom sets of spectra for F01, we see the spectra 
for i ma x = 17 is consistent until the octupole (Z=3) where it 
spikes. This suggest that due to the larger data size of F01 
it is more able to cope with the demands of increasing the 
series expansion. Therefore, at times throughout this sec- 
tion, we will focus our analysis on the RM map from the 
F01 catalogue with Z m ax = 16. 

In Figures 2, 3 and 4, we show the RM maps for the 
S81, B88 and F01 catalogues, respectively. We display only 
the /max=10 and 16 maps in order not to overload the 
reader with information. The r.m.s. values of 1Z for the S81, 
B88 and F01 RM maps (J max =16) arc 26.4, 23.5 and 21.5 
radm -2 , respectively. Maxima (large positive TV) are white, 
whereas, minima (large negative TV) are black. The two lim- 
iting scales enable us to see the progress of structure as 
the series expansion is extended to include higher modes. 
We can observe how feature at small I develop as the se- 
ries extends. Reassuringly, the main features in the ' m ax = 16 
maps are also present in the Z ma x=10 maps. The positions 
of the maxima and minima remain roughly unchanged. This 
is compelling evidence that the observed features are real. 
However, comparison of the maps from the three catalogues 
is inhibited by the temperature-colour scale varying from 
map to map. Therefore, for the Z ma x = 16 maps, we force the 
maximum and minimum scale to be \1Z\ = 100 rad m -2 ; the 
results of which are shown in Figure 5. 



( 




Figure 3. B88 catalogue with 644 sources. Top: Z m ax =10. Bot- 
tom: Z max =16. 




I 



\ * 



Figure 4. F01 catalogue with 744 sources. Top: (max =10. Bot- 
tom: Z max =16. 



From Figure 5, it is clear that the maxima at Z n ~270° 
and minima at l u ~ 90° are the dominant features in all three 
plots (here we use Z n to denote Galactic longitude to avoid 
confusion with the angular scale I). This maxima/minima 
pair corresponds to the large-scale magnetic field in the lo- 
cal Orion spur (sometimes referred to as an arm) . These two 
spots are displaced from the equator to negative Galactic co- 
ordinates. This asymmetry between the two hemispheres has 
been widely reported before (eg. Vallee & Kronberg 1975; 
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Figure 5. RM maps with identical temperature-colour scaling 
('max = 16). From top to bottom: S81, B88 and F01 catalogues. 

Frick et al. 2001); it is usually attributed to the local radio 
Loop I (the North Galactic spur). Such local distortions are 
associated with interstellar magnetised superbubbles with 
typical diameters of 200 pc (Vallee 1997). There is also a 
prominent maxima/minima pair towards the Galactic centre 
in the RM map formed from the F01 catalogue. The centres 
of the maxima and minima are at l u — 1° and I = 346°, 
respectively . On closer inspection, this feature is present in 
the RM maps from the other two catalogues. Finally, in the 
S81 RM map, there is a strong maxima at I ~ 50° in the 
northern hemisphere. This feature is only suggested in the 
other two maps. 

Cross-sections, along the Galactic equator, were taken 
of the RM maps (Z ma x=16) in order to further understand 
the features. These are shown in Figure 6. Ideally, with such 
a slice, maxima/minima locations should indicate the tan- 
gential direction to spiral arms and directional field changes 
should correspond to 71 = 0. However, local distortions and 
flaws in the map-making process, make this not entirely true. 
The orion spur location is clear for all three maps. Further- 
more, the maxima/minima pair towards the Galactic centre 
(described in the previous paragraph) is evident in all three 
cross-sections. However, the picture is hazy from Z n =30-50°: 
there is clear field reversal in S81 map; a hint of a reversal 



o 
o 




180° 90° 0° 270° 180° 

Figure 6. Galactic equator cross-section (imax — 

16). Dashed: S81 

catalogue; Dotted: B88 catalogue; and Solid line: F01 catalogue 
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Figure 7. F01 catalogue with both the dipole and quadrupolc 
removed (Z max = 16). 

in F01 map; and none in B88 map. The maxima and min- 
ima in the cross-sections could be attributed to the named 
inner arms (eg. the minima near the Galactic centre to the 
Norma arm), however, this seems quite speculative given the 
variation from catalogue to catalogue. 

It is clear that the Orion spur is the dominant feature. 
Since the associated maxima/minima pair is separated by 
180°, it will be the main source of the dipole (Z=l). More- 
over, we see from the spectra that the quadrupole (£=2) 
is also strong. Therefore, we remove both the dipole and 
quadruple from the RM map (Z ma x=16) compiled from the 
F01 data. The results of which are displayed in Figure 7. 
This enable us to view some of the smaller scale features. 
These details will be hard to explain solely from Galactic 
magnetic field models. It will be interesting to see, if these 
small-scale features persist with larger data sets. 

A study of the global Galactic magnetic field structure 
would benefit from the removal of local distortions. Conse- 
quently, we applied the method with the region containing 
Loop I removed (6 > 0°, < l u < 40°, 270 < F < 360°) 
(Ruzmaikin & Sokoloff 1977). However, the removed seg- 
ment was too large to successfully reconstruct the sky given 
the remaining sources. The lack of restrictions in the seg- 
ment meant large maxima/minima formed there. This high- 
lights one disadvantage of spherical harmonic analysis over 
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wavelet analysis that can be localised in both physical and 
wavenumber spaces. The removal of these local structures is 
useful for getting a clear picture of the Galactic magnetic 
field. In CMB foreground studies, however, these structures 
are essential components of a template. 

We now turn to the question of errors in the derived RM 
maps. In principle there are two distinct types of uncertainty 
that could arise in the analysis we described above. The first 
concerns : errors intrinsic to the measurement of 1Z and the 
second relates to errors resulting from sample selection. 

We have tackled the first type of error by removing 
sources with the most extreme values of 1Z (those with 
TZ > 300radm~ 2 ), on the grounds that these are least likely 
to be galactic. We expect the remaining experimental errors 
to be stochastic and therefore the process of extracting in- 
formation on large-scales (as we do) should be unaffected by 
any underlying noise. The spectrum of this noise should be 
flat and only dominate on scales where the "real" power is 
weakest. Experimental errors were addressed in some detail 
by Frick et al. (2001) in the construction of their catalogue 
and their versions of the two other catalogues used in our 
analysis. We feel it would be inappropriate to repeat such a 
detailed analysis here. 

The second type of error corresponds to the sampling 
uncertainties. This error could be computed reliably if we 
had a large number of independent samples. This is, of 
course, impossible but even in cases of non-repeatable ob- 
servations there are resampling techniques that can be used 
to make reasonable estimates of these errors. The purpose of 
resampling the data is to generate further sets with the same 
population distribution as the original. Small perturbations 
to the original data set will lead to this. For example, Ling, 
Frenk & Barrow (1986) apply a 'bootstrap' resampling tech- 
nique to estimate the sampling errors in the two-point cor- 
relation function estimated from galaxy and cluster redshift 
data. The bootstrap technique involves sampling N points 
(with replacement) from the original data set of N sources 
in order to create pseudo data sets. The variation over an 
ensemble of such resamplings is used to estimate the error 
in the statistic in question for the original data. We applied 
this method to the F01 catalogue with the intention of find- 
ing the sample errors in the RM map with Z max set to 16. 
In total, we produced 50 bootstrap samples from the origi- 
nal data set and from each of these constructed a RM map. 
The resampling technique is only designed to test the in- 
ternal variance of the true data set. Mean values obtained 
from the ensemble of pseudo data sets are not expected to 
be good estimators of the true mean values. Therefore, from 
the bootstrap samples, we calculate the standard deviation 
a at each pixel position p 



a(p) = 



\ 



E 



K^p) - (Tl(p)} 
49 



(18) 



where i corresponds to the bootstrap sample. From this, we 
constructed a signal-to-noise map where the signal is taken 
as 1 72. | and the noise is a. This map is displayed in Figure 
8. The map saturates at a signal-to-noise of unity so the 
regions with low signal to noise are more obvious. The mean 
signal-to-noise across the whole sky is 117. We found this 
value to be very stable as the number of bootstrap samples 
increased and this dictated the number of pseudo data sets 




Figure 8. Signal-to-noisc map constructed from 50 bootstrap 
samples of the F01 catalogue. Dark regions are where the esti- 
mated signal-to-noise exceeds unity. 



produced. Clearly, the majority of the sky is unaffected by 
sample errors and we feel reassured that the features seen 
are not the result of sample selection. It appears that the 
boundaries between positive and negative regions of the sky 
where the signal is correspondingly low are the most suscep- 
tible to sample selection. 

The spherical harmonic coefficients generated for 
all three catalogues with i max set to 16 are avail- 
able at http://www.nottingham.ac.uk/~ppxptd/rm_maps. 
Hopefully, this will enable our method to be compared with 
other techniques and allow the maps to used in the investiga- 
tion of observables affected by the Galactic magnetic field. 
Instructions on the generation of full-sky maps using the 
HEALPix package are also given at this address. 



5 CORRELATIONS WITH CMB MAPS 

As mentioned in the introduction, RM maps have a gen- 
eral importance beyond trying to map the Galactic mag- 
netic field. In what follows, we hope to display one particular 
function. In this section, we focus solely on the RM map pro- 
duced from the F01 catalogue with i max set to 16. Dineen 

6 Coles (2004) developed a diagnostic of foreground con- 
tamination in CMB maps. The method measured the cross- 
correlation between the RM of extragalactic sources and the 
observed microwave signals at the same angular position. 
In what follows, we seek correlations between the spheri- 
cal harmonic modes of the RM map and CMB-only maps. 
In doing so, we shall look at the phases of the (complex) 
coefficients of the modes from 1=2 to Z=16. Phase correla- 
tions have been used before to hunt for evidence of depar- 
tures in the CMB temperature field from a Gaussian random 
field (Coles et al. 2004). In Dineen, Rocha & Coles (2004) 
a certain form of phase correlation was found to be associ- 
ated with non-trivial topologies. Phase correlations between 
CMB and foreground maps have been sought before (Nasel- 
sky, Doroshkevich & Verkhodanov 2003; Chiang & Nasel- 
sky 2004), however, here we wish to emphasise the virtue of 
having independent probes of Galactic foreground contami- 
nation. 

In order to seek evidence of phase correlations between 
the RM map and CMB data, we turned to two WMAP- 
derived maps. Both were constructed in a manner that min- 
imises foreground contamination and detector noise, leaving 
a pure CMB signal. The ultimate goal of these maps is to 
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build an accurate image of the last scattering surface (LSS) 
that captures the detailed morphology. Following the release 
of the WMAP 1 yr data, the WMAP team (Bennett et al. 
2003), and Tegmark, de Oliveira-Costa & Hamilton (2003; 
TOH) have released CMB-only sky maps (see papers for de- 
tails) . We use the WMAP team's internal linear combination 
(ILC) map and the Wiener-filtered map of TOH. The latter 
was chosen since the map was found to be correlated with 
RM values in Dineen & Coles (2004). 

Two measures of phase association were used: the circu- 
lar cross- correlation coefficient R and Kuiper's statistic V . 
Both statistics will be evaluated at each scale 2 from 2 to 
16. If we let $rm and $cmb be the phases of the RM and 
CMB maps, respectively. Then, following Fisher (1993), R 
is defined as: 
i 

R(l) = r 1 COs($ m , RM - $m,CMB). (19) 

m— 1 

The expectation value of R is 0, and hence highly corre- 
lated phases are associated with large values of \R\. Kuiper's 
statistic is calculated from the available set of phase differ- 
ences ($ m ,RM — $ m ,cMB) at a given scale. First, the phase 
differences are sorted into ascending order, to give the set 
{0i, • • • , Op}- Each angle Qj is divided by 27r to give a set 
of variables Xj , where j = 1 ... p. From the set of Xj we 
derive two values Sp and Sp where 

S+=max|i-Xi,^-X 2 ,...,l-X p | (20) 
and 

5'-=max|xi,X 2 -i,...,X p -^^J. (21) 
Kuiper's statistic is then defined as 

V(l) = (5+ + Sp)- ^+0.155 + ^^. (22) 

The form of V is chosen so that it is approximately inde- 
pendent of sample size for large p. Anomalously large values 
of V indicate a distribution that is more clumped than a 
uniformly random distribution, while low values mean that 
angles are more regular. 

To access the significance of the values of 7? and V ob- 
tained from the comparison of the RM map with the two 
CMB maps, we make use of Monte Carlo (MC) skies with 
uniformly random phases. The statistics were calculated for 
10,000 MC skies contrasted with a further 10,000 MC skies. 
Thus, we are left with 10,000 values of R and V for each 
scale. 

The results from both CMB maps suggests that there is 
strong correlations between the phases at 2=11. For the ILC, 
the values of -R(ll) and V(ll) are greater than 99 percent 
of the MC values. Whereas, the values of -R(ll) and V(ll) 
corresponding to the TOH Wiener-filtered map are greater 
than 97 and 98 percent of the MC skies, respectively. In 
Figure 9 we plot the ILC map constructed with the ai tTn for 
2=11. Overlapping this image with that of the RM in Fig 
4, we can see that the central maxima/minima pair in the 
RM map are similar in location, size and shape to structure 
in the ILC image (but with colours reversed). This is prob- 
ably what determines the specific scale 2=11. Interestingly, 




Figure 9. Internal linear combination map constructed with a; j7n 
for £=11 only. 

2=11 corresponds to the scale that Naselsky, Doroshkevich 
& Verkhodanov (2003) found the greatest level of correla- 
tion between the ILC phases and those of the foreground 
maps. 



6 CONCLUSION 

In this paper, we have presented a new method to generate 
all-sky RM maps from uneven and sparsely populated data 
samples. The method calculates a set of functions orthonor- 
mal to the data set. With these basis functions, the spherical 
harmonic coefficients are calculated and converted into sky 
maps using the HEALPix package. The method was applied 
to three catalogues; S81, B88 and F01 catalogues. Maps from 
each catalogue showed evidence of the magnetic field in our 
local Orion spur, the North-South asymmetry attributed to 
radio Loop I and a maxima/minima pair close to the Galac- 
tic centre that possibly corresponds to the magnetic field of 
two inner Galactic arms. A RM map constructed from the 
S81 catalogue also had a prominent maxima at 2 n ~ 50° in 
the northern hemisphere. 

In Section 5, we showed the benefits a RM map has to 
CMB foreground analysis. Phase correlations were sought 
between RM maps and those of CMB-only maps derived 
from the WMAP data. For both the WMAP team's ILC 
map and the Wiener-filtered map of TOH, phases corre- 
sponding to 2=11 were found to be highly correlated. Nasel- 
sky, Doroshkevich & Verkhodanov (2003) found the same 
scale to display phase correlations when carrying out a sim- 
ilar analysis on the ILC map and foreground maps. Their 
detection of correlations at the same scale as our analysis, 
reaffirms that the RM catalogues are valuable independent 
tracers of CMB foregrounds (Dineen & Coles 2004). 

Modelling foregrounds will play a crucial role in CMB 
polarisation studies. Foreground contamination is expected 
to be more severe than in the temperature measurements 
(Kosowsky 1999). Consequently, superior templates for the 
individual foreground components are required. RM maps 
will help trace these components. Besides this, extrapolation 
of low frequency measurement of synchrotron polarisation 
to CMB frequencies has been shown to be complicated by 
Faraday rotation (de Oliveira-Costa et al. 2003). Again, this 
underlines the importance of developing templates of the 
Faraday rotation of the Galactic sky. Efforts to map the RM 
sky will be greatly enhanced by increased source catalogues 
for both extragalactic sources and pulsars within our Galaxy. 
This may enable the formation of a 3-dimensional image of 
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the Galactic magnetic field. Furthermore, attempts to map 
the RM sky will be enhanced by upcoming satellite CMB 
polarisation experiments which present unprecedented sky- 
coverage and resolution. 
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